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Abstract 

We include correlations between elemental abundances in a Monte Carlo 
statistical analysis of BBN predictions, which, along with updated reaction 
rates and an improved BBN code, lead to tightened constraints on Qb and 
N v . Observational upper limits on the primordial 4 He and D + 3 He fractions 
of 24% (by weight) and 10~ 4 respectively lead to the limits: 0.015 < Qb < 
0.070, and N v < 3.04. The former limit appears to be incompatible with 
purely baryonic galactic halo dark matter, while the latter puts qualitatively 
new constraints on neutrinos, and physics beyond the standard model. 
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The remarkable agreement, both qualitative and quantitative, of the pre- 
dicted primordial light element abundances and those inferred from present 
observations yields some of the strongest evidence in favor of a homogeneous 
FRW Big Bang cosmology. Because of this, significant efforts have taken 
place over 20 years to refine BBN predictions, and the observational con- 
straints they are related to. Several factors have contributed to the maturing 
of this field, including the incorporation of elements beyond 4 He in compari- 
son between theory and observation (i.e. |J), and more recently: an updated 
BBN code @, a more accurate measured neutron half life [Q] , new estimates 
of the actual primordial 4 He , D+ 3 He, and 7 Li abundances f|, ||, and finally 
the determination of BBN uncertainties via Monte Carlo analysis 0. All of 
these, when combined together ]7] yield a consistent picture of homogeneous 
BBN which is at the same time strongly constrained by observation. 

We have returned to re-analyze BBN constraints motivated by three fac- 
tors: new measurements of 

several BBN reactions, the development of an improved BBN code, and 
finally the realization that a correct statistical determination of BBN predic- 
tions should include correlations between the different elemental abundances. 
Each serves to further restrict the allowed range of the relevant cosmological 
observables flBaryon and N v . 

1. New BBN Reaction Rates: By far the most accurately measured 
BBN input parameter is the neutron half-life, which governs the strength of 
the weak interaction which interconverts neutrons and protons, and which 
in turn helps govern when this reaction drops out of equilibrium. Since this 
effectively determines the abundance of free neutrons at the onset of BBN, it 
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is crucial in determining the remnant abundance of 4 He . With the advent of 
neutron trapping, the uncertainty in the neutron half-life quickly dropped to 
less than 0.5% by 1990. Nevertheless, it is the uncertainty in this parameter 
that governs the uncertainty in the predicted 4 He abundance. The world 
average for the neutron half-life is now tn = 889 ± 2.1 sec ||, which has an 
uncertainty which is almost twice as small as that used in previous published 
BBN analyses j|, || [7j. We utilize the updated value in our analysis. 
Next, a recent measurement of 

7 Be + p — > 7 + 8 B suggests[|] a rate about 20% smaller at low energies 
than previous estimates. One might expect that at high values of 

7710 (defined by the relation Q B = .0036/i- 2 (T/2.726) 3 ^i x 10 10 , where T 
is the microwave background temperature today, and h defines the Hubble 
parameter by H = lOO/i km/ (Mpc sec)) lowering this rate would result in less 
7 Be destruction, which would increase the 7 Li abundance resulting from the 
decay of 7 Be after BBN. However, this is a subdominant destruction process 
for 7 Be. 

We find that the reducing the rate by 20% in our code alters the remnant 
7 Li abundance by less than one part in 10 5 ! 

Outside of these, we updated the Kawano code to use the reaction rates 
and uncertainties from in Smith et al [0. 

2. New BBN Monte Carlo: Because of the new importance of small 
corrections to the 4 He abundance when comparing BBN predictions and ob- 
servations, increased attention has been paid recently to effects which may 
alter this abundance at the 1% level or less. In our BBN Code several such ef- 
fects were incorporated, resulting in an r] 10 -independent correction of +.0006 
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to the lowest order value of Y p (the 4 He mass fraction). This is a change of 
+.0031 compared to the value used in previous published analyses j|, |5], [7J. 

This earlier value was based on correcting the lowest order value of Y p by 
an amount —.0025 0, based on the work of Dicus et al ||. The Dicus et al 
correction has two significant pieces: -.0013 from integrating the weak rates 
rather than using an expansion in powers of T to calculate \{n <-> p), and 

-.0009 from using the correct Coulomb correction, rather than simply 
scaling the neutron lifetime. This latter approximation 

incorrectly "Coulomb corrects" rates which do not feel the electromag- 
netic potential, such as ne + — ► pi>, and also ignores any temperature depen- 
dence. The remaining corrections — radiative, finite temperature, electron 
mass effects and neutrino heating — either 

effectively cancel (the first two) or are insignificant (the last two) ||. 

In the present code, more than half of the new correction 

is due to finer integration of the nuclear abundances. Making the time- 
step in the code short enough that different Runge-Kutta drivers result in 
the same number for the 4 He abundance produces a nearly 7710 independent 
change in Y p of +.0017 fK|. The other major change is the inclusion of 
M^ 1 effects [l]]. Seckel showed that the effects on the weak rates due to 
nucleon recoil, weak magnetism, thermal motion of the nucleon target and 
time dilation of the neutron lifetime combine to increase Y p by ~ .0012. 

The kinematics of nuclear recoil, are responsible for roughly 

25% of the increase which Seckel found [|1(J. Also included in the correc- 
tion is an small increase of .0002 in Y p from momentum dependent neutrino 



decoupling [O, |13 
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Finally, we have utilized a Monte Carlo procedure in order to incorporate 
existing uncertainties and determine confidence limits on parameters. Such 
a procedure was first carried out by Krauss and Romanelli |J, who chose 
the BBN reaction rates from a (temperature-independent) distribution based 
on then existing experimental uncertainties. Their procedure was further 
refined by Smith et al J/]], who both updated the experimental uncertainties, 
and utilized temperature dependent uncertainties in their analysis. Here 
we utilized the nuclear reaction rate uncertainties quoted by Smith et al 
(including the temperature dependent uncertainties for 3 He(a,7) 7 Be and 
3 H(a, 7) 7 Li) except for the updated reactions described earlier. Each reaction 
rate was determined using a Gaussian distributed random variable centered 
on unity, with a width based on the quoted 1 — a uncertainty from Smith et 
al. For the 

rates without temperature dependent uncertainties this number was used 
as a multiplier 

throughout the nuclear abundance integration. For the two rates with 
temperature dependent uncertainties the original uniformly distributed ran- 
dom number was saved and mapped into a new gaussian 

distribution with the appropriate width at each time step. 

While the Smith et al analysis cut off each distribution at ±2.6cr, we did 
not made such a restriction. Our code was designed to generate warnings, 
discard data, and reset random numbers if reaction rate values which were 
generated became nonsensical (i.e. negative). Warnings were generated only 
for the temperature dependent rates. For 3 H(a,7) 7 Li 1 warning per 4000 
BBN runs was generated, less than 1 warning per 30000 runs was generated 
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for 3 He(a, 7) 7 Be. 

The results of our updated BBN Monte Carlo analysis are displayed in 
figure 1, where the symmetric 95% confidence level predictions for each ele- 
mental abundance are plotted. Also shown are claimed upper limits for each 
of the light elements f|, |], based on 

observation^. This figure also allows one to assess the significance of the 
size of the corrections we have used in relation to the width of the 95% 
C.L. band for Y p , which turns out to be ~ .002. The total change in Y p of 
~ +.003 from previous BBN analyses conspires with the reduced uncertainty 
in the neutron lifetime, which narrows the uncertainty in Y p and also feeds 
into the uncertainties in the other light elements, to reduce the acceptable 
range where the predicted BBN abundances are consistent with the inferred 
primordial abundances. 

3. Statistical Correlations Between Predicted Abundances: While 
the introduction of a Monte Carlo procedure was a significant improvement 
in the evaluation of BBN uncertainties and predictions, the determination of 
limits on the allowed range of BBN parameters Qbaryon and N v based on com- 
parison of symmetric 95% confidence limits for single elemental abundances 
with observations, as has become the standard procedure, overestimates the 
allowed range. This is because the BBN reaction network ties together all 
reactions, so that the predicted elemental abundances are not statistically 
independent. In addition, the use of symmetric confidence limits is too con- 
servative. Addressing both of these factors is a central feature of our work. 
Figure 2 displays the locus of predicted values for the fractions Y p and 
2 Where the estimates differ, we have used the more conservative one. 
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D + 3 He/H for 1000 BBN models generated from the distributions described 
above for i] 1Q = 2.71 (figure a) and r) w = 3.08 (figure b). Also shown is 
the x 2 — 4 joint confidence level contour derived from this distribution, 
in a Gaussian approximation, calculating variances and covariances in the 
standard manner. The horizontal and vertical tangents to this contour cor- 
respond to the individual symmetric ±2er limits on Gaussianly distributed x 
and y variables. As can be seen, the distribution is close to Gaussian, but has 
deviations. Nevertheless, this approximation is useful to quantify the magni- 
tude of correlations and variances. We calculated the normalized covariance 
matrices at different values of 7710, and display the covariances in Table 1. As 
is evident from this table, as well as the figure, and as is also well known on 
the basis of analytical arguments, there is a strong anti-correlation between 
Y p and the remnant D + 3 He abundance. Thus, those models where 4 He is 
lower than the mean, and which therefore may be allowed by the existing 
quoted upper bound of 24% on Y p , will also generally produce a larger rem- 
nant D+ 3 He/H abundance, which can be in conflict with the quoted upper 
bound on this combination of 10~ 4 . This will have the effect of reducing the 
parameter space which is consistent with both limits, as we now describe. 

Because our Monte Carlo generates the actual distribution of abundances, 
Gaussian or not, we determine a 95% confidence limit on the allowed range 
of r] W (N p ) by requiring that at least 50 models out of 1000 lie within the 
joint range bounded by both the 4 He and D + 3 He upper limits, as shown in 
figure 2. This is to be compared with the procedure which one would follow 
without considering joint probability distributions. In this case, one would 
simply check whether 50 models lie either to the right of the 
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D + 3 He constraint for low 7710 (figure a), or below the 4 He constraint for 
high 7/10 (figure b). This is clearly a looser constraint than that obtained 
using the joint distribution. Finally, the procedure which has been used 
to-date, which is to check whether the symmetric 2cr confidence limit for a 
single elemental abundance crosses into the allowed region gives even a looser 
constraint, as can be seen in figure 2a. This, after all corresponds to checking 
whether only 25 models lie to the right of the D + 3 He constraint for low 7710 
(figure a), or below the 4 He constraint for high 7] w (figure b). 

In table 2 and figure 3 we display our results. Here we show the 95% confi- 
dence limits on r/i , as we have defined them above, and also using the looser 
procedures which ignore correlations. As can be seen, accounting for the 
correlations in the non-symmetric 95% confidence limit tightens constraints. 
Moreover, the impact of the procedure becomes stronger as the number of 
effective light neutrino species, N v is increased. Greater than 3.04 effective 
light neutrino types is ruled out only once correlations are taken into account. 

We also determined an upper limit on r] W using just 7 Li. Requiring 
7 Li/H < 2.3 x 1(T 10 §, [7J yields a limit r/ 10 < 5.27 . This is weaker than 
the 4 He limit, and there remains some debate about the actual observational 
upper limit on primordial 7 Li (i.e. changing 2.3 to 1.40 will lower the limit 
on 7/10 to 4.15.). Alternatively, we can use the bound on 7/10 derived above to 
set an allowed range of 9 x 10~ n — > 1.5 x 10~ 10 on the primordial values of 
7 Li, which should be compared with the observational estimates. 

4. Conclusions and Implications: The new 

constraints we have derived here on t/ 10 , and N u , taken at face value, have 
significant implications for cosmology, 
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dark matter, and particle physics. The limit on 7710 corresponds to the 
limit 0.015 < 0, Baryon < 0.070. (To derive this bound we required 0.4 < h < 
0.8, as is required by direct measurements and 

limits on the age of the universe.) Thus, if the quoted observational upper 
limits on the Y p and D+ 3 He/H are valid homogeneous BBN implies that: 

(a) The upper limit on flbaryon seems 

incompatible with all, or even galactic halo dark matter being purely 
baryonic. 

(b) The bound on the number of effective light degrees of freedom during 
nucleosynthesis is very severe, corresponding to less than 0.04 extra light 
neutrinos. This is a qualitatively different constraint than the previously 
quoted limit of 0.3 extra neutrinos. For example, it rules out any Dirac mass 
for a neutrino without some extension of the standard model because even a 
right handed component which freezes out at temperatures in excess of 300 



GeV will contribute in excess of 0.047 extra neutrinos during BBN i.e. ||14 
without extra particles 

introduced whose annihilation can further suppress its abundance. Even 
allowing 0.047 extra light neutrinos, the upper limit on a Dirac mass would be 



reduced to « 5 keV JL7|, [H] • Similarly, new light scalars are ruled out unless 
they decouple above 300GeV. A v T mass greater than 0.5 MeV with lifetime 
exceeding 1 sec. is also ruled out due to its effect on the expansion rate 



during BBN i.e. see [15], [16 1 - Also, neutrino interactions induced by extended 



technicolor at scales less than O(100) TeV are ruled out [ljj. Moreover, sterile 
right handed neutrinos |H| would be ruled out as warm dark matter as the 
lower limit on their mass would now be 0(lkeV). ( We will explore these 
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constraints further in subsequent work). 

Finally, having devoted considerable effort to accounting for the statis- 
tical uncertainties in BBN predictions, we must still stress that the largest, 
and most significant, uncertainties in the comparison of BBN predictions 
with observations come from the latter. Moreover, the uncertainties in these 
observational limits are dominated by systematic, and not statistical effects. 
Hence, the 95% confidence limits we derive must be qualified by the recogni- 
tion that their significance is really only as good as the observational limits 
are. Such limits cannot, at present, be taken to imply statistically invio- 
late constraints on neutrino parameters or Qb- In other words, systematic 
errors in the quoted upper limits on the inferred primordial light element 
abundances could allow the limits quoted here to be broadened. 

Nevertheless, the theory can be carefully tested. If, for example, bary- 
onic dark matter is found to make up the galactic halo or if neutrino mass 
measurements conflict with the above bounds, this would most likely imply 
that the quoted upper limits on 4 He , or on 3 He + D are flawed. This would 
be of great interest for stellar evolution studies. 

Indeed, the existing constraints from BBN are now so tight — requiring a 
primordial 4 He fraction in excess of 23.8% for consistency — that an agnostic 
view is prudent at the present time as to whether the constraints derived 
above will be satisfied or else whether observations will require revision in 
the inferred primordial abundance estimates. Finally, we note that inho- 
mogeneous BBN is not likely to alter this conclusion, as recent work has 



established. |2f |. 
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Table 1: Normalized Covariances 





Vio 


covariance of 


1.00 


2.00 


2.50 


3.00 


3.40 


4.00 


4 He vs D + 3 He 


-.71 


-.60 


-.49 


-.47 


-.44 


-.6 


4 He vs 7 Li 


-.23 


.04 


.12 


.25 


.30 


.25 



Table 2: Correlations & 7710 limits 



95% C.L. 
7710 range 


N v 


3.0 


3.025 


3.04 


3.05 


w/ corr. 


2.69 3.12 
2.65 <-> 3.14 
2.62 3.17 


2.75 <-> 2.98 
2.65 <-> 3.04 
2.63 <-> 3.10 


2.83 2.89 
2.69 2.99 
2.65 3.03 




2.69 <-> 2.95 
2.66 <-> 3.00 


w/out corr. 


sym. w/out corr. 
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Figure Captions 

Figure 1: BBN Monte Carlo predictions as a function of 7710- Shown are 
symmetric 95% confidence limits on each elemental abundance. Also shown 
are claimed upper limits inferred from observation. 

Figure 2: Monte Carlo BBN predictions for Y p vs D+ 3 He and allowed range 
for (a) r]i Q = 2.71, and (b) r) w = 3.08. 

In (a) a Gaussian contour with ±2a limits on each individual variable is 
also shown. 

Figure 3: Number of models (out of 1000 total models) 

which satisfy constraints Y p < 24% and D + 3 He/ H < 10~ 4 as a function 
of 7710, for 3.0,3.025,3.04,3.05 effective light neutrino species. Curves are 
smoothed splines fit to the data. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/94020 1 Ov 1 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/94020 1 Ov 1 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/94020 1 Ov 1 
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